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Abstract 

We calculate the c-axis kinetic energy difference between normal and super- 
conducting state for coherent and for incoherent interlayer coupling between 
Cu02 planes. For coherent coupling the ratio of the missing conductivity 
spectral weight to the superfluid density is equal to one and there is no vio- 
lation of the conventional sum rule, but for the incoherent case we find it is 
always greater than one whatever the nature of the impurity potential may 
be. To model more explicitly YE^CusC^-z around optimum doping, which 
is found to obey the sum rule, we consider a plane-chain model and show 
that the sum rule still applies. A violation of the sum rule of either sign is 
found even for coherent coupling when the in-plane density of electronic states 

depends on energy on a scale of the order of the gap. 
PACS numbers: 74.20.-z, 74.25. Gz 
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It has been proposed that the interlayer coupling along the c axis of a high tran- 
sition temperature (T c ) superconductor is incoherent, and the electronic kinetic energy 
along the c axis changes when the system enters the superconducting state. 0~@ Recently, 
Basov et a/.i have reported that there is a significant discrepancy between the superfluid 
density p s and the spectral weight missing from the real part of the c-axis conductivity 
N n — N s = 8 Jq+ duj &i c (uj) — G\ c (oj) , where u c is a cutoff frequency of the order of a band- 
width, in several high-T c cuprate superconductors such as optimally doped Tl 2 Ba2Cu6+ x 
(T12201). This implies that the conventional sum rule of Ferrel, Grover and Tinkham (FGT)i 
is violated. However, the spectral discrepancy becomes vanishing for YBa2Cu3 06.85 and dis- 
appears for the optimally doped YBa 2 Cu30 6 .9 5 (YBCO) crystal with T c ~ 93K as the c-axis 
response become more coherent with increasing oxygen content. Basov et al. also pointed 
out that there is no such discrepancy in the in-plane response for any cuprate. Moreover, 
for over-doped T12201, the sum rule discrepancy vanishes and a Drude-like peak develops 
in the conductivity for T > T c .@ Theses observations, therefore, suggest that for coherent 
interlayer coupling in the cuprate superconductors the conventional sum rule is obeyed. 

In this paper, we consider both coherent and incoherent c-axis coupling between Cu02 
planes. For the coherent case we find that the superfluid density remains equal to the missing 
optical spectral weight; in other words, it does not violate the FGT sum rule. The c-axis 
kinetic energies in the normal and superconducting state have the same value. For incoherent 
c-axis coupling the ratio of the missing area to the superfluid density is always larger than 
one in disagreement with some recent experiments. In YBCO the CuO chains play an 
important role in the electrodynamics and at the optimal doping a plane-chain model§0 
is needed to be complete. Here we use this model to investigate the c-axis conductivity 
sum rule. An algebraic calculation of the electronic kinetic energies is complicated and a 
numerical calculation is required although it can be reduced a lot in a special case, in which 
only the leading order in perturbation theory is kept. This case is particularly interesting 
because it has been shown to exhibit a pseudogap in the c-axis conductivity.0 Finally, we 
discuss the possibility that the FGT sum rule is violated in the plane-plane case when the 



in-plane density of states depends on energy even for coherent case. 

The Hamiltonian H for a cuprate superconductor with coherent c-axis coupling is H = 
H + H c , where H describes a (i-wave superconductor in a plane and H c = J2ia t±[ c iia c i2a + 
cf 2(T Ciia\ is a coherent interlayer coupling due to the overlap of electronic wave functions 
which is represented by t±; therefore, by coherent coupling we mean a tight binding like 
coupling along the c axis. It will not be necessary to treat t± as a constant in what follows. 
It can depend on an angle in the plane. For incoherent coupling the Hamiltonian is H' c = 
J2iaVi[cf la Ci2a + c&aCiio]) where Vi is an impurity scattering potential, so that impurity 
scattering mediates the c-axis hopping and an impurity average is implied.0Ei 

In the presence of an external vector potential A z , H c is modified to H C (A Z ) by the phase 
factor exp(— ieA z ) for c^c^o- and exp(ieA z ) for cf^Ciia- For the response to an external field, 
H C (A Z ) is expanded up to second order of A z to obtain the current j c = —5H C (A Z )/5A Z = 
jp + jd, where j p = —iedY,i a t±\cti a Ci2a — c^cua] and jd = e 2 d 2 H c A z with d the interlayer 
spacing. In linear response theory, (j c ) = [— U + e 2 d 2 (H c )]A z , where II is the current-current 
correlation function associated with j p and (H c ) is the perturbation of jd due to H c . The 
conductivity cx c (q, u) is given by 



a c q,u;j = - 



' H(^u)-e 2 d 2 {H c ) . (1) 



In the Matsubara formalism, 



n(q, u) = 2(ed) 2 T ^ ^ ^Tr [f G(k, u')f G(k, J + u)] , (2) 

u>> k 



and 



(H c ) =2TY,T. [^^(k, w)f 3 G(k, u)\ , (3) 

UJ k 

where f% is the Pauli matrix in the spin space, and G(k, u) is the Green's function in Nambu 
representation, namely, G(k, u) = — (iuiTo + ^kT"3 — AkTi)/(u; 2 + ^k + ^k) wiht ^k the in-plane 
energy and Ak the gap which has d x 2_ y 2 symmetry in the cuprates. 
The c-axis conductivity sum ru of the system is 



2 r°° 

- / dua lc (u) = -e 2 d 2 {H c ). (4) 

7T JO 

We use the unit such that H = c = k% = 1 and set the volume of the system to be unity. 
From the sum rule, the superfluid density p s can be written as 



fW-C r -i r -i 

jf du[al(u)-al c (u)]-47re 2 d 2 [(H c ) s -(H c ) n ], (5) 



where uj c is the cutoff frequency for interband transitions that H c does not account for. 

Since the difference between the superfluid density and the missing spectral weight is 
proportional to the kinetic energy difference between normal and superconducting state as 
seen in Eq. (|5]), it is necessary to calculate (H c ) s — (H c ) n to see if the FGT sum rule is 
violated by coherent c-axis coupling. For the normal state, 

(H c ) n = 4T££tiG (k,u;) 2 , (6) 

u k 

where G (\c, uj) is a normal state Green's function and t± may depend on k z and <fi = 
tan -1 (k y j k x ) . We assume a cylindrical Fermi surface with £ = k 2 /2m — /i, where /i is a 
chemical potential in the plane, and a <i-wave gap Ak = A(T) cos20fc. Then, we obtain 

K z 

where the integration range is limited by u c , and the density of states, iV(£), is approximated 
by a constant value N(0) around the Fermi energy. Later, we will discuss the effect of iV(£) on 
(H c ) and will note the possibility that the FGT sum rule may be violated even for coherent 
c-axis coupling. Since df(£)/d£ = —5(0 at zero temperature (T = 0), (H c ) n turns out 
to be — 4iV(0) J2k z I d<fi / (27r)t±. For a superconducting state with superconducting Green's 
functions G(k,u) and F(k,u), 



(H c y= AT^Y,^ [ G ( k > - F (K 



w k 



V 2 + 

The difference between {H c ) s and (H c ) n is of the order of (A(T)/c<j c ) 2 ; therefore, coherent 
c-axis coupling does not violate the FGT sum rule as long as u c » A(0) even if t± depends 
on cf). Note that the difference is largest at T = and vanishes as T — ► T c . 

The calculations for incoherent (impurity mediated) c-axis coupling proceed in the same 
way as before. Note that in this case j p = —iedj^ia Vi[ c iia c i2cT — cf 2r j c iio-] an d id = e 2 d 2 H' c A z , 
and an impurity configuration average is required. We derive the normalized missing spectral 
weight (N n — N s )/p s under assumption of a constant density of states and show that it is 
greater than one. 

The penetration depth A c can be calculated in two ways. Based on the Kramers-Kronig 
relation for the conductivity, we obtain A c , namely, l/47rA 2 = ]xm u) ^,o[utma c (0, ui)]. Alter- 
natively, using Eq.(H) we can also calculate A c (= l/ v /p7). Equate these two expressions of 
A c , then after integration over energy we arrive the formula as follows: 

(jV„ - N.) i i ^/#fc#,|y(&,6JI 2 [i - gggggg ] 

P. 2 + 2 E./dMWV(*»,«l' 7 ^ 7 ^ 3 j 

The second term in Eq. ([|) can easily be shown to be bigger than one half whatever the an- 
gular dependence of the impurity potential V(4>k, 4> p ) may be.il Thus the normalized missing 
spectral weight is always greater than one. In a simple model of impurity scattering,B0 for 
which \V((p k ,(f) p )\ 2 = \V \ 2 + |Vi| 2 cos20pCos20 fe , we found that (JV„ - N s )/p s > 1.58. This 
incoherent coupling model, therefore, does not agree with recent findings. The sum rule is 
one for YBCO around optimum doping indicating coherent c-axis coupling and less than 
one for the underdoped case. To treat YBCO around optimum doping more realistically we 
need to include the complications introduced by existence of the chains along the b axis. 

Penetration depth (X a a>)) experiments^ in YBCO have shown that both A a and A& are 
linear T at a low temperature and that a considerable amount of the condensate resides 
on the chains. To treat this case we need to consider a plane-chain coupling model.S0 We 
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assume the hybridization of Fermi surfaces between plane and chain arising through coherent 
coupling. For simplicity we also assume the gap in the chain has a <i-wave symmetry and its 
magnitude is of the order of that in the plane. The Hamiltonian for a coupled plane-chain 
system is H = £ k C£ fi, k Ck, where (7+ = (C^, Ci-kj., C£ kp C 2 -kM and 



/ 



fi k -A lk t(k z ) 
-Aik -6k -t(fc z ) 
t(k z ) Gk -A 2k 



V 



(10) 



/ 



-t{k s ) -A 2k -£ 2k 

where t(/c 2 ) = — to cos(k z d/2) for coherent coupling between plane and chain, £i( 2 ) is the 
energy dispersion in the plane (chain), and A^ 2 ) is a gap of the plane (chain). We point out 
here that the conclusion we make later does not depends on the simple form of t(k z ), and 
that Ai(2) and £i( 2 ) depend only on k x and k y . 

The Hamiltonian of the plane-chain coupling model is also decomposed into two parts, 
H = Hq + H c . Hq is for the superconductivity in the plane-chain coupling system and 



its eigenvalues can be reduced to ±E± = ±\/e± + A k , where e± are normal state energy 



dispersions e± = (^ + ^ 2 )/2 ± y(^i — ^ 2 ) 2 /4 + t{k z ) 2 with A lk = A 2k = A k for simplicity. 
Extensive work on this Hamiltonian can be found in Ref.i 

In order to calculate the linear response of the system to the external electromagnetic 
field, we modify H c with the phase factor mentioned before and follow the same procedure 
to derive the current j c = j p + jd- Then, we obtain 



Jp 



edto 



J2sin(k z d/2)C+a 1 ®t C k , 

k 



(11) 



where d/2 is the distance between a plane and a chain and a% is a Pauli matrix in the 
plane-chain space, and jd = e 2 (d/2) 2 H c A z , where H c = J2kt{^z)C^d- 1 <g> f 3 (7 k . The c-axis 
conductivity for q = 0, <r c (0, u>), of the sysyem is also derived to be cr c (0, cu) = (i/oj) n(0, uj) — 



e\d/2f{H c ) 



where 



n(0,cj) = (edt /2) 2 Tj2J2 sin2 ( k z d / 2 ) TT [ & i ® TbG(k, o;')o-i ®f G{k,u' + u) 

u 1 k 



;i2) 
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and 

(H c ) = ^T^^cos 2 (A;,d/2)Tr[ ( j 1 ® f 3 6(k,u;)£i ® f 3 G(k, cj)] , (13) 

w k 

with G(k,r) = -(T[£ k (r)C£(0)]), which is a (4 x 4) matrix. The Green's function 6(k,w) 
is given by G(k,u) = (iu — h^)^ 1 . We emphasize that the Hamiltonian in this model is 
qnite different from the nsnal macroscopic tunneling Hamiltonian,0 for which, for example, 
Gi 3 (k, r) = — (^"[C , ikT( r )C2kt(0)]) i s n °t allowed because each layer is independent (as is 
the case for the previous coupling model); however, it is possible in the present model 
because of the hybridization through the chain between the two Fermi surfaces of plane and 
chain. Introducing a unitary matrix U which diagonalizes hk, one can show Gy(k, o>) = 



YL=i UimU+j/ (iu - E m ), where E m = ±E± if A ik = A 



2k- 



(H c ) becomes complicated and the energy dispersion in the chain is quite different from 
that in the plane so that a numerical calculation is required to see if the difference between 
(H c ) s and {H c ) n is negligible. However, since to in Eq. ( |T3"D is assumed small we may expand 
(H c ) in terms of to and keep only the leading order, which is £q. This case includes only 
interband Hamiltonian but is still very interesting as it can exhibit a c-axis pseudogap.0 In 
this approximation with A lk = A 2k = A k and for special case, {H c ) s becomes 

^' -^^wrirm^kwy (14) 



Note that (H c ) s in Eq. ([14]) is almost same as (H c ) s in Eq. (H) for the simple coherent 
coupling case except that now ^ and d/2 appears rather than d in t(k z ). One, therefore, 
may expect that 5{H C ) will vanish to order (A(T)/u; c ) 2 . It is obvious that S(H C ) is identically 
zero along the nodal lines, and 5{H C ) is largest along the anti-nodal directions. Since £i = 
k 2 /2m — n and £ 2 = k 2 y j2ra — /i, we introduce £ = £i, then £ 2 = £ sin(0) 2 — /icos(0) 2 . If 
= 7r/2, then £ 2 — 6; therefore, it can be seen that 5 (-^c) 0=71-/2 is of the order of (A/u c ) 2 . 
For = 0, £ 2 = —\l and it can be shown that 



H c )l =0 ~ 2N(0)J2tl r rf£[-^tanh(A ; 

^ J-w c l t, 2 - /r 2T 



tanh( — ) . (15) 



(e 2 -/i 2 )VF + ^ V 2T 
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Now, the leading order of 5(H C )^ =Q changes to (A/yu) 2 . It is possible to show that for an 
arbitrary (ft, as long as fi and oj c » A, S(H C )^, is negligible, and consequently, the FGT sum 
rule is not violated in the plane-chain coupling model. 

In a numerical calculation for the general case without the above simplification, we have 
computed 5{H C ) / {H c ) s , which is the fractional change in kinetic energy. We have taken 
£2 = ky/2m — // for the chain energy dispersion, where // is a chemical potential in the 
chain. For simplicity, we also assume // = fi. It has been shown that // — //(<< y) may 
correspond to the pseudogap seen in the c-axis response of over- doped YBCO;0 however, 
it makes no difference in the numerical evaluation of 5(H C ) / (H c ) s . We choose T = 12K, 
A(T) = 20meV, t = 2meV, /i = 500 meV and uj c = 400meV. We found that S(H C )/(H C ) S 
becomes more negligible as we increase the summation range of the Matsubara frequency to. 
For \uj\ < 200ttT, 5(H c )/(H c ) s ~ 8.8 xlO" 3 , and for \u\ < 2000vrT, 5(H C )/(H C ) S ~ 5.4 xlO" 3 . 

One may consider a plane-plane coupling through a chain. In order to investigate the 
c-axis kinetic energy for such a coupling, one needs to replace £2 and A 2 with £x and Ai, 
respectively. For the hopping amplitude, t(k z ) can be simply changed to t(k z ) 2 because 
the plane-chain and chain-plane distances are the same and equal to d/2. Then, one can 
algebraically show that S(H C ) is as negligible as before. It is also possible to see that S(H C ) 
has a symmetry with respect to £1 <-> £2 and Ai <-> A2; in other words, 5{H C ) for the chain- 
plane coupling is the same as that for the plane-chain coupling. Therefore, it implies that 
8(H C ) along the c axis is conserved for coherent coupling. 

So far we have taken the density of states as a constant: iV(£) = iV(0) (— uo c < £ < u c ) 
and concluded that the difference of the c-axis electronic kinetic energies between normal 
and supercondcuting state is negligible. Now we would like to consider the effect of iV(£) 
on the sum rule when it is a function of £ to illustrate possible changes. If it varies strongly 
with £, it clearly cannot be approximated by N(0). We taylor-expand N(£) up to £ 2 near 
£ = 0. Then, jV(0 = N(0) + £0iV(OMIo + (£/2) 2 <9 2 iV(0/<9£ 2 |o. At T = 0, (H c ) n of 
Eq. (|?p does not change; however, Eq. (g) for (H c ) s changes due to (£/2) 2 iV" (0), where 
N" (0) = d 2 N(£)/d!; 2 \o. Assuming t± in Eqs. (0) and (H) does not depend on 0, we obtain 
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5(H c )/(H c ) n ~ (8iV(0))- 1 iV"(0)A(0) 2 ln(^/A(0)). (16) 

Note that this correction can have either sign depending on the sign of the second derivative. 
For JV"(0)/iV(0) ~ l/cu 2 c , 5(H c )/(H c ) n ~ or 1 - 65 , where x = w c /A(0), because ln(x)/x 2 ~ 
3.-1.66 when x » L if N"(0)/N(0) ~ l/(A(0)w c ), then 5(H c )/(H c ) n ~ x" - 65 ; thus, <5(# c ) 
is considerable. For this to be the case AT' (0) needs to exhibit variation on an energy scale 
of order A rather than u c . In a realistic model a Taylor expansion about £ = may not be 
accurate but our calculations serve to illustrate the main point. Violation of the FGT sum 
rule of either sign can result from an energy dependence in the in-plane electronic density 
of states iV(£). The exact amount depends on details and cannot be known without a 
specific knowledge of the band structure involved. In-plane dynamics gets reflected in c-axis 
properties. 

For coherent interlayer coupling between Cu0 2 planes the superfluid density is equal to 
the missing optical weight; the FGT sum rule is satisfied. This applies even in more realistic 
model for YBCO around optimum doping such as the plane-chain model with two atoms per 
unit cell. On the other hand incoherent c-axis coupling mediated through impurity scattering 
gives a sum rule which is always larger than one and is in disagreement with experiment. To 
get the sum rule to be less than one as is observed in underdoped YBCO and other systems 
such as optimally doped T12201, it may be necessary to go to more exotic non- Fermi liquid 
pseudogap model for the in-plane motion as discussed recently by Ioffe and Millis.0 Their 
arguments however do not apply to optimally doped T12201 because this system does not 
show a pseudogap. Their pseudogap argument that leads to the cancellation of G(k, u) 
and Go(k, u) contribution to the ratio of missing area to superfluid density making it one 
half instead of one for the preformed pair model was made for coherent c-axis coupling, but 
we find it also applies to the incoherent case. El Another interesting model for the in-plane 
dynamics is the "mode" model of Norman et a/.0 introduced from consideration of ARPES 
data. In more conventional models a sum rule violation of either sign can also be obtained 
if there is a strong energy dependence to the density of states near the Fermi surface on the 
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scale of a few times the gap. 
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